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Abstract TUM 


Starting from the monopolar representation of the Rao-Wilton-Glisson (RWG) 
functions, an approximate inverse (AI) of the RWG Gram matrix is constructed, where 
only 3x3 matrices for all individual triangles have to be inverted. The 


inversion of the Gram matrix is, however, only approximate, since it results as the 
inverse of a product of singular matrices, which is not easily obtained. The AI is shown 
to be an effective preconditioner for the RWG Gram matrix and is in particular useful 
for ill-conditioned Gram matrics on distorted meshes. 


Index Terms: Gram matrix, Rao-Wilton-Glisson functions, preconditioner, matrix 
inversion 
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Motivation 


Gram Matrix Inversion: Part of Many Integral Equations 
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Discretization of Calderón projector and preconditioner 
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Testing function mappings from Buffa-Christiansen to Rao-Wilton-Glisson (RWG) 
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marching-on-in-time scheme for the time-domain 
combined field integral equation, TAP 2013. 
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Accurate MFIE discretization scheme with RWG functions 
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Kornprobst and Eibert, “An accurate low-order discretization 


scheme for the identity operator in the magnetic field 
and combined field integral equations,” TAP 2018. 
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In this presentation: focus on RWG Gram matrix 


Outline 


The Magnetic Field Integral Equation 
Improved MFIE: Weak-Form Identity Operator Representation 


Approximate Inverse Construction of the RWG Gram Matrix 


Numerical Results 


The Scattering Problem — Magnetic Field Integral Equation (MFIE) TLITI 
Huygens' Principle for a Perfect Electrically Conducting (PEC) Object 


PEC boundary condition 

nxe=0, n-h=0 oF We 
Superposition of incident and scattered fields * $ à So 
nxe"*-nxe*?-0 


y iot PEC object 
n x h" +n x h? =n x h"? 


nxe=0 


The Scattering Problem — Magnetic Field Integral Equation (MFIE) TLITI 
Huygens' Principle for a Perfect Electrically Conducting (PEC) Object 


PEC boundary condition 


nxe=0, n-h=0 p we 
n 
Superposition of incident and scattered fields X $ [ So 
inc sca _ tot (ee ? 
nxh'+nxh*=nxh \ zero total field ( 
Equivalent free-space problem N inside the X 
N . \ 
> Electric surface current densities ją = n x h'^' (Love currents) d MERE d 
E 5 
oll AG: = = Js 
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» Total field h'^' in terms of j, 


j; =n x h" +n x hä? 


The Scattering Problem — Magnetic Field Integral Equation (MFIE) TLITI 
Huygens' Principle for a Perfect Electrically Conducting (PEC) Object 


PEC boundary condition 


nxe-90, n:h-0 hinc hic 
n 
Superposition of incident and scattered fields 1 $ [ So 
inc sca _ tot (eee P 
nxh™+nxW@=nxh \ zero total field ( 
Equivalent free-space problem: h**? in terms of js N inside the N 
JA source region — 
e j 
inc lif e S a e E 
nxh'* z [17 - X] j, p 


Tic Hs, Js(r’) &(r, r^) ds’ 
Kj,:=nx fs, VGo(r, r’) X js(r’) ds’ 


MFIE Discretization 


Classical Discretization with RWG functions 


> Div-conforming RWG basis functions f 


> electric current densities j, ~ ¥ Pli]; in vector i 
» Div-conforming testing functions f not in the dual space of n x h 
> Conditioning mostly influenced by the Gram matrix Ggg 
> Poor accuracy due to wrong testing of 7 and 'K operators 
» Tested incident magnetic field vector h with the entries 


[A]m = nno : n x h"*(r) da 


» Tested 'K operator with the matrix entries 


[Kee] mu Jf B.C) “nx [[ VGo(r, r’) X By (r’) da’ da 


MFIE Discretization 


Classical Discretization with RWG functions 
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Div-conforming RWG basis functions fj 


> electric current densities j, ~ ¥ Pli]; in vector i 
Div-conforming testing functions f not in the dual space of n x h 
Conditioning mostly influenced by the Gram matrix Ggg 

Poor accuracy due to wrong testing of J and K operators 

Tested incident magnetic field vector h with the entries 


[Alm = nno : n x h"*(r) da 


Tested J operator with the matrix entries with strongly singular Dirac delta kernel 


[Gsp] am^ Jf B.C) . Af Bur )8(r, r’) da’ da 


MFIE Accuracy Problems 
A Long-Known Problem 
> Inaccuracies of Classical MFIE due to wrong testing functions 


> Most severe for sharp-edge objects 
> here a small wedge 
> bistatic radar crosssection o shows —9 GB rel. error for MFIE 


> for comparison: EFIE has —22 dB error on same mesh 
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Intermezzo — 90? Rotation 


Weak-Form Rotation of RWG-Discretized Currents 


> Div-conforming RWG basis $ for js 


> Rotated RWG functions a =n x B 


> not feasible, line charges etc. 


> Rotation expanded and tested by RWG functions 


Mapping of n x f functions 


on f functions 











MFIE with Weak-Form Identity Operator Representation TUM 


Two Subsequent Weak-Form Rotations 


> Identity operator is replaced by two 90° weak-form rotations and one 180° rotation 





> Weighted addition with standard Gram matrix; improved MFIE: 








MFIE with Weak-Form Identity Operator Representation TUM 


Two Subsequent Weak-Form Rotations 


> Identity operator is replaced by two 90° weak-form rotations and one 180° rotation 





> Weighted addition with standard Gram matrix; improved MFIE 


> Inverse Gpp required, obtained iteratively 


> diagonally-preconditioned conjugate-gradient solver? 
> should be solved more accurately than the scattering problem itself 


> increased condition number for deformed meshes 


MFIE Accuracy Problems — Solved! 


A Long-Known Problem 


» Improved identity operator discretization labeled as new MFIE 


» Same example as before: a wedge 
> bistatic radar crosssection c shows —9 dB error for MFIE 
» for comparison: EFIE has —22 dB error on same mesh 
» improved MFIE solution resembles EFIE and has —22 dB error 
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MFIE Accuracy Problems — Solved! TUM 


A Long-Known Problem 


> Improved identity operator discretization labeled as new MFIE 


> Same example as before: a wedge 


[1Ggg — 1 Gpo G5; po * Kgg] i — ^ 


> bistatic radar crosssection o shows —9 dB error for MFIE 
» for comparison: EFIE has —22 dB error on same mesh 
» improved MFIE solution resembles EFIE and has —22 dB error 





» This motivates us to study the Gram matrix inverse! 
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Computation of RWG Gram Matrix 


General Structure, Sparsity 


» Calculation of RWG Gram matrix entries 


[Gee] mn 7 // Br) : Bur) da 


> 5N entries in Gram matrix, sparse! 





» example wedge with N — 1284 and 3? angle 




















Computation of RWG Gram Matrix 
Triangle-Wise Integration 1 
» Calculation of RWG Gram matrix entries 


[Gee] „n = ff B) - Bur) da 


» nth pair of triangles is split into triangles n1 and n2 


» analytical integration, no quadrature 


[Gpp] mn 7 I] Bor) d Bu(r) da + I Bor) : p«r) da 


> N RWG functions B — 2N monopolar RWG functions B 


-plg.. 
Geg =P GggP 


Computation of RWG Gram Matrix 
Triangle-Wise Integration 2 


T 
» RWG Gram matrix via monopolar representation Gee =P Gij P 
> Ggg is NXN, Gjj is 2N X 2N 
» mapping matrix P goes from each edge to two triangles; must be tall 
> each row has only one entry, a factor +1 for the correct edge-triangle index 
> back-mapping from triangles to edges via PT 


> P and P! is just a notation for the summation of matrix entries 


» Block structure of monopolar RWG Gram matrix, diagonal elements G; are 3 x 3 


W 
W 
I 
. (~ 
o 
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Inverse of RWG Gram Matrix 


Inverse of a Matrix Product? 

» RWG Gram matrix via monopolar representation 
—-plg.. 

Gog =P Gjj P 


» Inverse of multiplication of (full-rank) matrices A, B 


(AB)! - BA"! 


Inverse of RWG Gram Matrix 


Inverse of a Matrix Product — or the Pseudo-Inverse? 
» RWG Gram matrix via monopolar representation 


-plgG.. 
G5; - P Gj;P 


» Inverse of multiplication of (full-rank) matrices A, B 
(AB)! - B-1A7! 
> Obviously, P is tall and thus singular — inverse does not exist 


» Try the pseudo-inverses instead? 


(AB)! - B'A' 


Inverse of RWG Gram Matrix 


Inverse of a Matrix Product — or the Pseudo-Inverse? 
» RWG Gram matrix via monopolar representation 


-plgG.. 
G5; - P Gj;P 


» Inverse of multiplication of (full-rank) matrices A, B 
(AB) - B-1A7! 
» Obviously, P is tall and thus singular — inverse does not exist 
» Try the pseudo-inverses instead? 
(AB)! - B'A' 
» Only valid under some circumstance which are neither fulfilled (for Moore-Penrose) 


» Try it anyhow, see what happens 


Approximated Inverse of RWG Gram Matrix 


Approximation, Pseudo-Inverse Not Working 
» Pseudo-inverses of mapping matrices are easy 
p! - 1pT pi! -ip 
2 i 2 
» Inverse of monopolar RWG Gram matrix calculated block-wise, G; is3 x3 
G! N 0 


-1 _ . 
Cig=|0 o c 


> Approximate inverse (AI) of Geg 
AI 1pT frz-1 


PE BB 


» Al has same structure as Ggg 


Comparison of Approximated Inverse to Gram Matrix Inverse TUM 


Obviously an Approximation 


Structure of Gee Structure of G5; with rel. threshold 107? 


(sparse, 5N entries) (fully populated) 
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Al as Preconditioner — Numerical Results 
Singular Value Spectrum 


» Sphere with 999 unknowns, singular value spectrum of 


> standard cond(Ggg) = 3.13 
> diagonally-preconditioned cond(G 5°) =2.5 


> Al-preconditioned cond(G55) = 1.28: lowest condition number 
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Al as Preconditioner — Numerical Results TUT 


Iterative Solver Convergence 

> hundreds of random right-hand sides b considered 
> conjugate gradient solver 

> solution error: ||r||z = ||x°® — x*||2 (CG solution x*8, real solution x*) 


> Improved solution accuracy, faster convergence 





1284 unknowns wedge 999 unknowns sphere 
stopping criterion 107!° stopping criterion 10? 
Matrix = its. avg. |[rll2 Matrix = its. avg. |Irll2 
Gpp 33 57.10? Gpp 12 18.1077 
6 27 58-10 Go 10 47.1974 
Gj5Ggg 16 14-10? GjjGgg 6 25:107 


Conclusion 


V 


Inverse of a Gram inverse important for many integral equations 


» example for improved MFIE discretization with RWG functions was presented 


V 


Analysis of RWG Gram matrix setup as multiplication of three matrices 


V 


Direct inverse not that simple due to pseudo-inverse properties 


» Approximate inverse as preconditioner increases accuracy and halves CG iterations. 
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